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Abstract We study the inhomogeneous generalization of a l-dimensional AKLT spin 
chain model. Spins at each lattice site could be different. Under certain conditions, the 
ground state of this AKLT model is unique and is described by the Valence-Bond-Solid 
(VBS) state. We calculate the density matrix of a contiguous block of bulk spins in 
this ground state. The density matrix is independent of spins outside the block. It is 
_i diagonalized and shown to be a projector onto a subspace. We prove that for large 

t^ , block the density matrix behaves as the identity in the subspace. The von Neumann 

r^ ' entropy coincides with Renyi entropy and is equal to the saturated value. 



Keywords AKLT ■ Density Matrix ■ Entanglement • Valence Bond Solid 
PACS 75.10.Pq ■ 03.65.Ud • 03.67.Mn ■ 03.67.-a 

1 Introduction 

Quantum entanglement is a fundamental measure of how much quantum effects we 
can observe and use to control one quantum system by another, and it is the primary 
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resource in quantum computation and quantum information processing [9], [49]. The 
entanglement of quantum states, particularly related with spin systems has been at- 
tracting a great deal of interest, see for example [5], [S], [18], [24], [26], [28], [32], [37j . 
[43] . [46] , [53] , [58] , [64] . [65] . [70] . Quantum entanglement can be quantified by the 
von Neumann entropy and Renyi entropy of a subsystem, as discussed in |18) . [21j . 
[22] , [36] , |38] , [66j , |71) . An area law of the von Neumann entropy in harmonic lattice 
systems has been proposed and studied in [13], [14] . [58] . Entanglement properties also 
play an important role in condensed matter physics, such as phase transitions [55], [56] 
and macroscopic properties of solids [2S], [62) . Extensive research has been undertaken 
to understand quantum entanglement for correlated electrons, interacting bosons as 
well as various other systems, see [3], [7], [10], [11], [H], [S], [H], [19], [27], [33], [34], 
[35], [39], [10], [42], [11], [45], [H], [51], [52], [5l], [57], [59], [60], [61], [62], [67], [68], 
[69], [72] for reviews and references. 

In this paper we study the inhomogeneous generalization of the AKLT spin chain 
model. The original homogeneous model was introduced by Affleck, Kennedy, Lieb and 
Tasaki [T], [5] with all bulk spins being the same. The inhomogeneous generalization 
was introduced in [41] in which spins at different lattice sites may take different values. 
The conditions of the uniqueness of the ground state were discussed and correlation 
functions in the ground state were obtained. The unique ground state is known as 
the Valence-Bond-Solid (VBS) state [6], [41], which is very important in condensed 
matter physics. The VBS state is closely related to Laughlin ansatz 47, and fractional 
quantum Hall effect 6_. It enables us to understand ground state properties of anti- 
ferromagnetic system with a Haldane gap [30]. Universal quantum computation based 
on VBS states has also been proposed [63| . 

The density matrix of a contiguous block of bulk spins (we call it the density matrix 
later for short) of the homogeneous AKLT model has been studied extensively in [18] . 
[23] ■ [38] ■ [41] ■ [64] ■ [71| . It contains information of all correlation functions [S], [37] . 
[38], [71] • Moreover, the density matrix was shown [18], [38] to be independent of the 
size of the chain and the location of the block relative to the ends. It was diagonalized 
in [71] and proved to be a projector onto a subspace. 

It was conjectured in [71) that the structure and properties of the density matrix 
is generalizable to: (i) the inhomogeneous chain; (ii) higher dimensional lattices; (iii) 
general graphs. In this paper, we shall prove the first part of the conjecture, i.e. that 
for the 1-dimensional inhomogeneous AKLT model. The general version of the inho- 
mogeneous model was first studied in [41]. In order to write down the Hamiltonian 
and conditions for the uniqueness of the ground state, we first introduce notations. 
Consider a system of a linear chain of A'^ bulk spins and two ending spins. By Sj we 
denote the vector spin at site j with spin value Sj . Then we associate a positive integer 
number to each bond of the lattice and denote by My (My = Mji) the bond number 
between sites i and j. They must be related to bulk spins by the following relation 

2Sj=M,„i,,+M,,,+i (1) 

with 250 = Mqi and 2S'/vf+i = M/^t^at _|_i for ending spins. The condition for solvability 
of relation ([l]) is 

N+l 

Y.{~iys,^o. (2) 



Solution to relation ([1} under condition ([2| is 

3 

M,.,j+i^2Y^{-iy-'Si>l. (3) 

1=0 

More details can be found in 41,. 

Now we defined the Hamiltonian of the inhomogeneous AKLT model as 

N Sj+Sj + i 

H = Y. Y. C7j(j,j + l)Pj(j,j + l). (4) 

Here the projector Pj{j,j + 1) describes interactions between neighboring spins j and 
j' + l, which projects the bond spin Jjj + i = Sj + Sj+i onto the subspace with total spin 
J {J = Sj + Sj^i + 1 — Mjj^i , . . . , Sj+ Sj^i) . An explicit expression of Pj (j, j + 1) was 
given in [4T]. Coefficient Cj{j,j + 1) can take arbitrary positive value. This Hamiltonian 
(|4]) has a unique ground state (VBS state). 

If we consider a block of L contiguous bulk spins as a subsystem, then it is possible 
to pick up the corresponding terms from Q which describe interactions within the 
block. The collection of these terms is referred to as the block Hamiltonian: 
k+L-2 Sj + Sj+i 

Hb= Y. Y. C7jO-i + l)Pj(j,i + l). (5) 

j=k J=Sj + Sj + i + l-Afj,j + i 
Here the block starts from site k and ends at site k + L — 1. The degeneracy of ground 
states of the block Hamiltonian ([5} is {AI^._i j, + l)(A/^,_|^^_x k+L + !)■ 

The density matrix p of the whole chain is a projector onto the unique VBS ground 
state (see (O). Because of entanglement with spins outside the block, the density matrix 
Pi of the block as a subsystem (see (|17p ) will no longer be a pure state density matrix 
as p is in general. We shall prove that the density matrix p^ is a projector onto a 
(Mj;_i fc + l)(Mj,_|_i_ijj_|_/, + l)-dimensional subspace of the complete Hilbert space 
associated with the block (see Sections[5]and[7]l . It turns out that the block Hamiltonian 
Hfi defines the density matrix p^ completely in the large block limit L ^ oo. The zero- 
energy ground states of the block Hamiltonian Hf, span the subspace that the density 
matrix p^ projects onto. So that p^ can be represented as the zero-temperature limit 
of the canonical ensemble density matrix defined by Hf,: 

PL = lim f — ^TTT' L ^ CO. (6) 

In the zero-temperature limit, contributions from excited states of Hi, all vanish and 
the right hand side of (|6} turns into a projector onto the ground states of the block 
Hamiltonian. 

As main subjects of the paper, we will construct eigenvectors and derive expressions 
for corresponding eigenvalues of the density matrix. We will show that the density 
matrix is a projector. The paper is divided into three parts: 

1. We calculate the density matrix using the Schwinger representation (Sections (2] [3]). 

2. We construct eigenvectors, derive eigenvalues and prove that the density matrix is 
a projector. (Sections |4l O [S]! 

3. We investigate the structure of the density matrix in the large block limit. As 
characteristic functions of quantum entanglement, the von Neumann entropy and 
the Renyi entropy are obtained in the limit (Section [7|. 



2 Unique Ground State of the Hamiltonian 

We start with the ground state of the Hamiltonian Q. It is given in the Schwinger 
representation by the VBS state [6], [41] 

N M- ■ 1 

|VBS)^[](a]&]+i~6]at^,) "'^' jvac), (7) 

where a* , b* are bosonic creation operators and |vac) is destroyed by any of the an- 
nihilation operators a, b. These operators satisfy [aj,aj] = [&i,fej] = Sij with all other 
commutators vanishing. The spin operators are represented as 57^ = alb-j, S'7 ~ bla-j, 

Sj = ^{a-aj — bjbj). To reproduce the dimension of the spin-Sj Hilbert space at site 
j, an additional constraint on the total boson occupation number is required, namely 
^(ajttj + bjbj) = Sj. More details and properties of the VBS state in the Schwinger 
representation can be found in [6], [8], [41] , The pure state density matrix of the VBS 
ground state is 

^ |VBS){VBS| 
^ (VBS|VBS) • ^ '' 

Normalization of the VBS state is 

N+l 

n (2^^- + !)• 

(VBSIVBS) = -jp^ . (9) 

See [71] for more details. 



3 Density Matrix of a Block of Bulk Spins 

We take a block of L contiguous bulk spins as a subsystem. In order to calculate 
the density matrix of the block, it is convenient to introduce a spin coherent state 
representation. We introduce spinor coordinates 

(m,«) = [cos-e*^,sin-e~*'^V < e < tt, < (^ < 27r. (10) 

Then for a point f2 = (sin6'cos(j!), sin^sin(;/>,cos6l) on the unit sphere, the spin-S co- 
herent state is defined as 

(ua* + vb^] 
\n) = ^ ' ' |va.c). (11) 

Here we have fixed the overall phase (a f/(l) gauge degree of freedom) since it has no 
physical content. Note that pi|) is covariant under SU{2) transforms (see Section[71). 
The set of coherent states is complete (but not orthogonal) such that [4], [23) 

5 
'^^±1 I dn\Q){Q\= J2 \S,m){S,m\^l2s+i, (12) 



where \S,m) denote the eigenstate of S and Sz, and /2S+1 is the identity of the 
(25' + l)-diniensional Hilbert space for spin-S'. 

Now we calculate the density matrix of a block of L contiguous bulk spins in the 
VBS state ((Tj. By definition, this is achieved by taking the pure state density matrix 
([8]) and tracing out all spin degrees of freedom outside the block: 



PL = rro4,...^fc-l,fc+L,...,iV,iV+l [p] 



1 < fc, fc + -L - 1 < iV. 



(13) 



Here the block of length L starts from site k and ends at site k + L — 1. pi^ is no longer 
a pure state density matrix because of entanglement of the block with the environment 
(sites outside the block of the spin chain). 

Using the coherent state representation (|lip and completeness relation (|12p . p^ 
can be written as |38], [71] 



PL 



(14) 



fe-l N+l 

n n '^^^ 

j=0 j=k+L 



fc-2 N 



n n [^(i-^.-^i; 



j=0 j=k+L 



Mi 



B^YBS l) (VBS l\B 



k+L-l 

n 

3=k 



{2Sj + 1) 
An 



N+l 



N 



1 — f2j ■ j/y-|_i) 



M, 



j=o 



Here the boundary operator B and block VBS state IVBS/^) are defined as 



B = {uk-ibk - Wfe-iflfc) ■'= i'*" {ak+L-l^k+L - h+L-l-u-k+L) 



M, 



k + L-l,k + L 



k+L-2 



AT- 

|VBSi)= II {a]b]^,-b}a]^i) "^' |vac) 



(15) 
(16) 



j=k 



respectively. Note that both B and |VBS/,) are SU{2) covariant (see Section[7]). After 



performing integrals over S7j (j — 0,1, . 



,k~2,k+L+l, 



,N, N+l) in the numerator 



and all integrals in the denominator, the density matrix p^ turns out to be independent 
of spins outside the block. This property has been proved for homogeneous AKLT 
models in |18) . [38) . [71) . Therefore, we can re-label spins within the block for notational 
convenience. Let k = 1 and the density matrix takes the form 



with 



PL = 



[](M„-+i + l) 
j=o 1 



1[{2S, + 1)\ 



, d^od^L+i-B^|VBSL)(VBSi|B 



i=i 



B' = [uob[ - voa[j (^fflL^L+l - o^wl+iJ 



Ml 



L-1 



ivBs^>=n(">]+i-^]4+i 



Mi, 



I vac). 



(17) 



(18) 
(19) 



The remaining two integrals of (|17|l can be performed, but we keep its present form 
for later use. 



4 Ground States of the Block Hamiltonian 



The block Hamiltonian ([5]) with the re-labeling k = 1 reads 



L — l Sj + Sj + i 

Hb = Yl Y. Cj{j,j + l)PjiJ,j + l). (20) 



Now we define a set of operators covariant under SU{2): 

4 = (^«I + "&!) ''^'' {4b[ - &I4) '^~' ("4 + <) 



where J_ = ^(Mqi - Ml^l+i), J+ = ^{Mqi + Ml,l+i) and \J-\ < J < J+. These 
operators (|2ip act on the direct product of Hilbert spaces associated with site 1 and 
site L. Then the set of ground states of (|2Up can be chosen as 

|G;J,^>=4|VBSi>, J=\J-\,...,J+. (22) 

Any state |G; J, fi) of this set depending on a discrete parameter J and a continuous 
unit vector Q \s a. zero-energy ground state of (|20[) . To prove this we need only to 
verify for any site j and bond {j,j -\- 1): (i) the total power of a- and h- is 25^-, so 
that we have spin-Sj at site j; (ii) — ^(Mj_ij- -I- Mj+i.j+2) < •^j.j+l = '^'1 + ^j+l — 
2(Afj"-lj+Mj_|_i ,j_|_2) by a binomial expansion, so that the maximum value of the bond 
spin Jjj+i is ^(Afj-ij -I- Afj+ij-|_2) — Sj + Sj+i — Mjj^i (from SU{2) invariance, 
see [6]). Therefore, the state |G; J, /?) defined in (|22|) has spin-Sj at site j and no 
projection onto the Jjj+i > Sj + Sj+i — Mjj^i subspace for any bond. 

The set of states {|G; J, /?)} with the same J value are not orthogonal. We also 
introduce alternatively an orthogonal basis in description of the degenerate zero-energy 
ground states of (|20|) . For notational convenience, we define 

X,,, . ^^:^i=;_, 4 - ("^)^-^"'(''¥-" (23) 

^(J + Af)!(J-M)! •^" y/{S + my.{S-m)\ 



These two variables transform conjugately with respect to one another under SU{2) 
(see Section (Tjl. Xjj\/ has the following orthogonality relation 

^ 47r 

■00 is a spin state creation operator such that 



V'Ll^ac) = |S,m). (25) 



The operator Aj defined in (|2ip can be expanded as (see [31], [71] ) 

^t ^ j (^-^+^+^)K^-+j)_K^+-^)K-^-+^)' x: x.M (26) 

mi-i-mL—M 
?7ii,?7ii:. 



where (^Moi,mi; ^Mi_i_|-i,7n2^| J, M) are the Clebsch-Gordan coefhcients. Note that 
tb\ . , and i/)I . , are defined in the Hilbert spaces of spins at site 1 and 

site L, respectively. We realize that the particular form of the sum over mi and tti^ in 
(|26|) can be identified as a single spin state creation operator 

mi+mL=M 

(27) 

This operator !^l > r acts on the direct product of two Hilbert spaces of spins at site 1 
and site L. It has the property that 

I'^mI''^'')^ ® |vac)L = I J, M)i^L. (28) 

If we define a set of degenerate VBS states {|VBSi(J, M))} such that 

|VBSL(J,Af)>=i^t^,|VBSL), J=|J-|,...,J+, M = -J,. ..,.], (29) 



then these (Mqi + 1)(A^l,L+i + 1) states (|29p are not only linearly independent but 
also mutually orthogonal (see Appendix [A} ■ Furthermore, any ground state |G; J, /?) 
can be written as a linear superposition over these degenerate VBS states as 



|G; J, n) 



(J+ + J + 1)!(J- + J)\{J+ - Jy.{-J- + J) 

2J + 1 
,/ 
• ^ Xjm\VBSl{J,M)), (30) 

M=-J 

and vice versa. More details can be found in [31) . [71) . 

Therefore, as seen from (|30|l . the rank of set of states {|G; J, /?)} with the same J 
value is 2 J + 1 and the total number of linearly independent states of the set { | G ; J, /?) } 
is ^ ,^| , |(2J + 1) = {Mqi + l)(Mi^_|_i + 1), which is exactly the degeneracy of the 

ground states of (|20p . So that {|G; J, /?)} forms a complete set of zero-energy ground 
states. The orthogonal set {|VBSi(J, M))} differs from {|G; J, O)} by a change of basis, 
so that it also forms a complete set of zero-energy ground states. 

5 Diagonalization of the Density Matrix 

It was shown in [71] for the homogeneous AKLT model that the density matrix is a 
projector onto a subspace. This statement is also true for the inhomogeneous model. 

Now we propose a theorem on the eigenvectors of the density matrix ([17p . The ex- 
plicit construction of eigenvectors yields a direct diagonalization of the density matrix. 
The set of eigenvectors also spans the subspace that the density matrix projects onto. 

Theorem 1 Eigenvectors of the density matrix pi^ | (J7| ) with non-zero eigenvalues are 
given by the set {\G;J,Q)} (Ml . or, eqmvalently, by the set {\VBSl{J,M))} [Ml . 
i.e. they are zero-energy ground states of the block Hamiltonian Hi, II20[) . 



We prove the theorem by re-writing the density matrix p^ (|17l) as a projector in 
diagonal form onto the orthogonal degenerate VBS states {\VBSl(J,M))} introduced 
in p9l). 

Take expression (|17|l and integrate over Qf) and Ol+i using binomial expansions 
and r\ dx(l + a;)"(l - a;)" = m!»! 2'"+"+i. Then we have 

L-1 

''' = — — ^ ^ U Jl ^ 

n(2S,+l)! ^=° "^^^ 

(VBSi|(6L)*^^-^+^-«(ai)«(ai)^'^«^-P(6if. (31) 

The particular combinations of bosonic operators appeared in pi|) are recognized up 
to a constant as spin creation operators "0 , . , , , , and ■(/',., , , , 

(see (|23[l ) at site 1 and site L, respectively. They commute with all bond operators 

f i i + + \^h.3 + l 

I tslb,- I ]^ — '^7"'^ 7+1 ) ' ®° that we could simplify the right hand side of (ISlfl using 

definition (|27|l and the following identity: 

l,L{vac|VlJ^,^^^ ij^^^^_p®ViA/x,.i,+i,g-iA/z„x + i 

Moi ^'^L.L^l 

11 11 

= E E \2^'^oi,i^Moi-p)i{i^Moi,-M^i^p\ 

1111 

(S>\^ML^L+l,q- -ML^L+l)L{-ML^L+l,q- -Ml^L+i\ 
J+ J 

= Y. Y. \j^m)^^l{j,m\ 

J=\J-\ M=-J 
J+ J 

= E E *'jMl^ac>i,i(vac|*'jM. (32) 

J=\J-\ M=-J 

The resultant final form of density matrix p^ is 
L-1 

PL^^ Moi!Mi,i+i! J2 E '^jA/lVBSi){VBSil<2/jM 

[](2S, + 1)! J=|./-|M=-J 



L-1 



= ^4 Moi\Ml,l+i^. J2 J2 |VBSi(J,M))(VBSi(J,M)|. 

(33) 

The set of degenerate VBS states {|VBSi( J, M))} with J = | J_ |, . . . , J+ and M = 
— J, . . . ,J forms an orthogonal basis (see Appendix [X)) . These (Mgi + 1)(M/, ^,+1 + 1) 
states also forms a complete set of zero-energy ground states of the block Hamiltonian 
(|20|) . So that in expression p3p we have re- written the density matrix as a projector 
in diagonal form over an orthogonal basis. Each degenerate VBS state |VBSi(J, M)) 
is an eigenvector of the density matrix, so as any of the state |G; J, Q) (because of the 
degeneracy of corresponding eigenvalues of the density matrix, see Section IS}. Thus we 
have proved Theorem 1. 



6 Eigenvalues of the Density Matrix 



Given the diagonalized form (|33|l . eigenvalues of the density matrix p^ can be derived 
from normalization of degenerate VBS states. We obtain an explicit expression for 
eigenvalues in terms of Wigner 3j-symbols in this section. 

Using expansion (|3Up and orthogonality (|24p . normalization of jVBSi(J, M)) can 
be calculated from integrating the inner product of jG; J, Q) with itself over the unit 
vector Q such that 



— / dniG; J, Q\G; J, Q) 
An J 

{J+ + J + ly.jj- + j)\{j+ - J)\{-J- + J)! 

(2J + 1)! 



(34) 



(VBSi(J,M)|VBSi(J,M)). 



This expression p4p also states that normalization of the degenerate VBS state is 
independent of /? and/or M (see Appendix 1X1 for another proof). 

Let's consider the integral involved in (|34|l . Using coherent state representation 
(|11|) and completeness relation (|12p as before, we obtain 



— / dn{G;J,n\G;J,n) 



(35) 



1 

4-K 



n 



{2Sj + 1)! 



47r 



dn 



L 



L-l 



m¥ 



1 — f2j • fij^i ) 



Mi 



-(i-ni-OL) 



J+-J V 



-[l + Q^.Q) 



\M, 



J- + .J 



[V 



I + O-Ql) 



-J-+J 



Now we expand [^(1 — i^i ■ ^j)] '^ in terms of spherical harmonics as in [23] . [38], 

[m 



-(i^a,-n,) 



Mij 



An 



Mi. 



M,j + 



-^A(z,Af„) J2 YUn,)YiUn^) 



(36) 



;=o 



10 



with 



X{l,M,, 



(-l)'My!(My+l)! 



(37) 



(M,j ~ iy.{M,j + 1 + ly: 

Then integrate over H and from ^2 to /^l-Ii the right hand side of (|35p is equal to 

L 

47rJ|(2Sj + l)! 
J=l 



L-1 

j=i 



{J- + J+ 1)( J+ - J + 1)(- J- + J + 1) 



L-l 



A/< J+-J J< ; i„ ifj 

EEEE E E 



■ \{la,J+ - J)A(/^, J- + J)\{lp, -J- +J) dQi / dfiL 

•V'^™(^l)il„,mJ^l)V]„™,(^l)ii:„,(^L)V5:,„.„(^L)y4,,„,(^L). (38) 

Where we have M< = mm{Mjj+i, j = 1, . . . , L — 1} and J< = min{J^ + J, ~J- + J}, 
both being the minimum of the corresponding set. Now we carry out remaining integrals 
in p8)) using 

dnYi^„,{f2)Yi^^„rAO)Y,^„rgm 



{21 + 1) {21a, +1) {21 + 1) lUo^lAl I la. l0 



(39) 



The result after integration can be further simplified by applying the following orthog- 
onality relation 



V(2;,+i)( ' ^" ^Mr '" '' 

^ — ' ^ \m ma rap I \mm, 






(40) 



where ( " I , etc. are the Wigner 37'-symbols. 

ym rua m/s J 

So that finally expression (|38} is equal to 



Y[{2s, + iy. 



L-l 
M< J+-J ,/< 

EEE 



(J- + J + 1){J+ - J + 1){-J- + J + 1) 



L-l 



nA(;,A'/,j+i) 



J=l 



A(/a, J+ - J)A(/^, J_ + J)X{lp, -J- + J) 



{2l + l){2l^ + l){2lp + l)fl^Y^ 



(41) 
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The explicit value of ( " '^ ) is given by 



I la Ip 





(42) 



(-1) 



g l{2g~2l)\{2g-2la)\{2g-2lpy- 



(25 + 1)! 



(g-l)\{g^la)\{g~lp)V 



li I -{-la + 1(3 = 2g (g £ N), otherwise zero. 

Combining results of (|33}, (|34|) and (|4H) . we arrive at the following theorem on 
eigenvalues; 



Theorem 2 Eigenvalues A{J) (J = \J-\, . . . , J+) of the density matrix are indepen- 
dent of n and/ or M in defining eigenvectors (see 1122]) and li29\) ). An explicit expression 
is given by the following triple sum 



A(J) 



L-l 

3 = 1 

L 

[](2S, + 1)! 

.7=1 



Moi!Mi,i+i!( VBSl{J, M)\VBSl{J, M)) 



(43) 



(2J+l)!Moi!Mi,i+i! 



{J+ + J+ 1)!(J_ + J + 1)!(J+ - J + 1)!(- J- + J + 1) 



M< 


J+-.J J< 


> 


>; >; 


I 


i„=0 i^=0 



L-l 



nA(i,M: 



j,i+iJ 



j=i 



A(Za, J+ - J)A(/^, J_ + J)\{ll3, ~J- + J) 



{2l + l){2la + l){2lp + l)(lYS) 



We shall emphasize at this point that given eigenvalues (|43p . both von Neumann 
entropy 



J+ 



Sv.N = -Tr[pi,\npi,] = - ^ (2J + l)yl(J)lnyl(J) 

J=\J-\ 



(44) 



and Renyi entropy 



Sn=Y^^n{Tr[pt]} = ^ln\ ^(2J + 1)^"(J) 



(45) 



J=\J-\ 



can be expressed directly. 
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7 Density Matrix in the Large Block Limit 



Consider the homogeneous AKLT model with spin-S at each site. In the limit L -^ oo 
(that is when the size of the block becomes large) , we learned from [16] , [29) , [38] , [66] 
that the von Neumann entropy reaches the saturated value Si, jv = In (S + 1) . It was 
also proved in [71] that the density matrix (denoted by poo in the limit) takes the form 

1 



Poo 



(5 + 1) 



iHs+ir 



><Po 



(46) 



where I/g^iy is the identity of dimension [S + 1) and •Poo is an infinite dimensional 
matrix with only zero entries. In this section, we generalize these properties in the 
limiting case to the inhomogeneous model. 

Let's apply the density matrix p^ p7[) to the state |G; J, i7) (|22p and get 



PL\G;J,n) 

L 

j=o 1^ 

L (47r)2 

3=1 



(47) 



df2odnL+iB''\VBSL){VBSL\BA''j\VBSL)- 



Using the coherent state representation (|lip and completeness relation (|12p , the factor 
{VBSl\BA^j\VBSl) in ((47| can be re- written as 



{VBSl\BAUyBSl} 



(48) 



(2gj- + 1) 
47r 



J ifldoA n [li^-^r^: 



n 

3 = 1 
■ (UQVl — VQUl) "^[UUl+VVlj {UlV]^—Vl 



Mi 



U*l) 



3 + 1 J 
J+-J 



We plug the expression (|48|l into (|47|l . Using transformation properties under 5*17(2) 
and binomial expansion, the integral over Qq yields that 



dOa \u},b\ - vla\\ {uavi - vqui) ° 



4tt 



Moi + 



- (uial +vibl 



t^ 



Moi 



J 



(49) 



Similarly we can perform the integral over 0[^^i. Then using expansion (|36|l and 
orthogonality of spherical harmonics, other integrals over fij with j = 2, . . . , L — 1 in 
can be performed. As a result, the following expression is obtained from (|47|l : 



M< 



1=0 



L-1 

J]a(/,m„-+i; 



KUf2)\VBSL). 



(50) 



13 



The operator K](0) involved in (|50p is defined as 



Moi 



kI{Q)= j dQidOLPli^l- ^L)iuia\+vih\] °\uu*i+vv*i)'''^'^ (51) 

/ * * * *\J+--J I * , *\-J-+J f t , ,t A*''^^-^^^ 

It is expressed as an integral depending on the order / of the Legendre polynomial 

Piini-OL)- 

There was no ambiguity in defining the large block limit in the homogeneous AKLT 
model [71) . However, in the inhomogeneous model we must specify the behavior of 
ending spins in the large block limit. So we define the large block limit as when L — > oo, 
the two ending spins approach definite values, namely, Mgi — > 5*-. and M^^^+i — * S+. 
Then we realize from (|37|) that as L — > oo, J^,_i X{l,Mjj^i) -^ 5iq. Therefore only 
the first term with Z = is left in (I5U|I . So that we need only to calculate the limiting 

Kim 

/c 
dQidQi^ (uior^ + fi&J ) (uui + vvij ' 

■ {ulvl - vlul) + (uul+vvl) " ■ iuLa\ + VLb\\ .(52) 

Here both J_ and ,7-f take the limiting values ^{S-—S+) and -^iS^+S^), respectively. 
It is useful to know transformation properties of the integrand in (|52|) under SU{2). 
The pair of variables (u, v) defined in (|10|) and bosonic annihilation operators (a, b) in 
the Schwinger representation both transform as spinors under SU{2). That is to say, 
if we take an arbitrary element D € SU{2) (2x2 matrix), then {u,v), etc. transform 
according to 

':)-K")- '=" 

On the other hand, {u*,v*), {—v,u), (a',fe^) and (~6, a) transform conjugately to 
{u,v). That is to say {u*,v*), etc. transform according to 

::)-D. (■■:), ,«, 

The integrand appeared in Kl{f2) (|52p . as well as Aj in (|21[) . boundary operator B' 
in (|18[) . etc. all transform covariantly under SU(2), i.e. those expressions keep their 
form in the new (transformed) coordinates. 

These transformation properties (|53p . (|54p can be used to simplify the Kl{f2) 
integral (see [71] for more details). We first make a SU{2) transform 



under the part of the integral H52p over fii . Then after integrating over Oi by a binomial 
expansion, we make an inverse transform in using D^^^ — DJ,j^ , consequently (|52p is 
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put in a form with a single integral over /?/, remaining: 

Klin) = ^i^ {ua\ + vblY'^' (56) 

/ t * t *\'^^^'^ r * *\-J-+J/ t t\'^+ 

Similarly, we make another SU(2) transform using 

It 

and integrate over i?/,. At last we make an inverse transform using D„ = D„, the 
final form is 

^o(^)- (5_AT(l + l) 4- (58) 

This expression states that {|G; J, /?)} is a set of eigenvectors of the density matrix as 
L —* oo. Let's denote the density matrix in the limit by poo- Then (|58p leads to the 
result (see ([50]l ') 

P^|G;J,^)=(^-^^i^^^|G;J,^). (59) 

We find from ((59} that the limiting eigenvalue yloo = ,o ^^iWc ^_i-i is independent 
of J. Any vector of the (S_ + 1)(S'+ + l)-dimensional subspace spanned by the set 
{|G; J, fi)} is an eigenvector of poo with the same eigenvalue ,„ +i'i/g +i'i • Therefore 
Poo acts on this subspace as (proportional to) the identity J(s_+i)(S++i)- So that we 
have proved explicitly that the density matrix takes the form 

P°^ ^ (S- + l)(5+ + l)-^(^-+i)(^++i) ® '^°°' ^^°^ 

which is a generalization of (|46p in the large block limit. In addition, we also derive from 
the eigenvalues that the von Neumann entropy •S'lj.Af = — X/ /— I / I (^"^ + l)^oo In Aoo 

coincides with the Renyi entropy 5'^' — y^ l'^ i X/ /— l / [(^^^ + 1)^Sd f a-^id is equal 
to the saturated value In [(5*- + 1)(S'+ + 1)]. 

8 Conclusion 

We have studied the density matrix p^ of a block of L contiguous bulk spins as a 
subsystem in the inhomogeneous AKLT model. The Hamiltonian is given by Q , which 
has a unique ground state under condition ([2}. The ground state is described by the 
VBS state ((7} in the Schwinger representation. The density matrix pj^ (|17[) of the 
block is obtained by taking trace (|13p of all spin degrees of freedom outside the block. 
The structure of the density matrix has been investigated both for finite and infinite 
blocks. It has been shown that the structure and properties of the density matrix are 
all generalizable to the inhomogeneous model, as conjectured in j71j . 
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For finite block, two mathematically rigorous results have been established as The- 
orem 1 and 2. In Theorem 1 we constructed eigenvectors of the density matrix with 
non-zero eigenvalues. These eigenvectors \G;J,n) defined in (|22p . or |VBSi(J, M)) 
defined in (|29p equivalently, are proved to be the (Mqi + 1){M]^ j^^i + 1) zero-energy 
ground states of the block Hamiltonian (I20|l . Using orthogonal basis {|VBSi(J, M))}, 
in Theorem 2 an explicit expression (|43|l for corresponding eigenvalues in terms 
of Wigner 3j-symbols is derived. Non-zero eigenvalues A{J) with J = |J_|,|J_| -|- 
1, . . . , J+ depend only on J and is independent of O and/or M in defining eigenvec- 
tors. The density matrix (|33p is a projector onto the subspace of dimension (Mgi -|- 
l)(A^L,L+l + 1) spanned by the set of eigenvectors {|G; J, O)} and/or {|VBSl( J, M))}. 

In the large block limit L ^> ex, ending spins approach definite values, namely, 
Mqi -^ S- and Mj^j^j^i -^ S+. All non-zero eigenvalues Aoo become the same (|59p . 
The infinite dimensional density matrix poo (|60p is a projector onto a (S_ -I- 1)(S'+ -I- 1)- 
dimensional subspace in which it is proportional to the identity. The von Neumann 
entropy S^jv coincides with the Renyi entropy S^ and is equal to the saturated value 
In [{S- + 1){S+ + 1)]. In the limit the Renyi entropy is a independent, which behaves 
quite differently from the XY model where the Renyi entropy has an essential singu- 
larity as a function of a (see [36) , [21j . [22]). 
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A Orthogonality of Degenerate VBS States 

The set of degenerate VBS states {|VBSi,(J, Af)), J = | J_|, ..., J+, M = -J, ..., J} introduced 
in I I29I I are mutually orthogonal. To show this, it is convenient to introduce the total spin 
operators of the subsystem: 



ij ij ij 

Stot=J2''l''^' ^tot = Y.''l''^' SU=J2^{ala,-blb,). (61) 

j=i j=i j=i 

First we show that the set of operators {•S'lti , ST ^ , S? ^ } commute with the product of valence 
bonds, i.e. 

L-l L-l 

These commutation relations 1621 1 can be shown in similar ways. Take the commutator with 
Srti first. We re- write the commutator as 

L-l 



L 
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L-l 



J = l 






(63) 



Then using commutators [ai,al] = Sij and [bi,fe.] = Sij, we find that 
[S+ + 5+, , . (atfot _ fet^t )Mj.,j+i] 

= aU-Mi i+i)a''..,(a''.bt., - b\at,A^'3,J+i-'^ 
W, ,M,- ,+ia't(a'':b't _ fot t jA/j^j+i-i 

j + l J.Ji-i jV J j + 1 ] J + 1' 

= 0. (64) 

Therefore [Stot'Il^Ii ("-Wj+l " '']4+l)*^'''"^^l = "• ^'^ ''^H' ^^^ 1^^™ "^^^ 

[bi,(a]6]+i - ftjaj+j'^^^-^+i] = -Mjj+ia]^^ia]b]^^ - fe]a]_^ J^-f^.^+i"!. In a parallel way, 

we find that the commutator with SIZ^ also vanishes. Next we consider the commutator with 

■^tof 



lStovl[^-]bl^-b]al^r 



■ 3 J + J- J J + J-' ' 



L-1 



(albl - blal)^'^^ ■ ■ ■ [5| + 5|+i, (a]b]+, - b]a]+ j'^^-.+i] 



j=i 



(65) 



In the right hand side of II65I I. the commutator involved also vanishes because 

= -\a\ai - b^'.b, + a'^..,ai+i - b'^..,b,+i,iah'^.., - fe'ta''^ , , )*^J.3+il 

2 3-' 3 ■' J + 1 -'+ J + 1 J + -^' ^ J J + 1 J J + 1' J 

+ a]^,[a,+^,{a]b]^, - b]a]^,)'^h.3 + ^\ - b]^,[b,+^,{,a]b]^, - b]a]^,)'^h.3 + ^\ 

= (66) 

Substituting ^^ into JSU, we obtain [S^^of 11^^/ ("j^J+i ~*]"]+i)^^^'^+^] = ^- ^°^ ™^ s'^'^" 
show that the state |VBSi^(J, M)) is an eigenstate of Sf + and the square of the total spin 

^tot = \ ('5'tof^t'^t + ■^tof^ttit) + (■S'tot)^ ""^^^ eigenvalues M and J{J + 1), respectively. Using 
the commutation relations 11621 1. we can show that 

L-i 



S±JVBSi(J,M)> = J|(a>J+i - ;>]at_^i)"^-^ + U5f +S±)|J,M)i,i,|vac>2,...,i_i 
Sf„JVBS£(J,M)> = HH^Vi - bJat^i)*^-^ + USf +S£)lJ,M>i,ilvac>2,...,i_i. 



(67) 
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Then from the definition of the state |VBSl(i/, M)) and the following relations: 



(5+ + 5+)|J,M>i,i = ^{J T M){J ±M + 1)1 J, M ± 1), 

{SJ= + S£)|J,M>i,i) = Af|J,M)i,i,, (68) 



we obtain 



S±JVBS£(J,M)) = ^{J :f M)(J ±M + 1)|VBS£(J, M ± 1)), 
■St^„JVBSz,(J, Af)) = Af |VBSi,(J, M)> (69) 

and hence S2^^|VBSi(J, A/)> = J(J + l)|VBSi(J, M)>. It is now proved that |VBSl(J, A/)) 
is an eigenstate of S'^^^ and S^^^ with eigenvalues M and J(J + 1), respectively. Therefore the 
states with different eigenvalues (J, Af) are orthogonal to each other. 

With the introduction of total spin operators of the block Kil l, it is straightforward to 
prove that normalization of the degenerate VBS state |VBSl(J, A/)) depends only on J and 
is independent of M . We prove the statement as follows: 



(VBSi,(J, M ± 1)|VBSl(J, M ± 1)> 



(VBSl(J, M)|STt5± JVBSz,(J, M)) 
1 



(J T M){J ± Af ~^^^^^^'^' *^)l(^tot - (-^tot)' ^ 5t^ot)|VBSi(J, M)> 
= (VBSi(J,Af)|VBSL(J,A/)>. (70) 

Here we have used the fact that |VBSl(J, Af)) is the eigenstate of S? ^ and Sf ^ with eigen- 
values J{J + 1) and Af , respectively. 
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